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1.  Statement  of  the  Problem  Studied 

The  major  objectives  of  this  project  are:  (1)  to  develop  a  generalized  parameterization 
method  in  three  dimensions  for  the  generalized  inverse  problem,  (2)  to  build  the 
theoretical  basis  for  model  structure  identifiability,  and  (3)  to  develop  a  new 
methodology  for  observation  network  design  that  unifies  model  structure  complexity, 
parameter  identifiability  and  model  application  reliability. 

2.  Summary  of  the  most  important  results 

We  have  developed  two  global-local  optimization  methods  for  three-dimensional 
parameter  structure  identification  in  groundwater  modeling.  Parameter  structure 
identification  is  formulated  in  terms  of  solving  a  generalized  inverse  problem  (GIP), 
which  allows  for  a  determination  of  an  appropriate  level  of  parameter  structure 
complexity,  and  the  identification  of  its  pattern  and  the  associated  parameter  values. 
The  structure  complexity  is  determined  by  calculating  the  parameter  structure  error 
measured  in  the  prediction  space  or  management  space,  while  the  structure  pattern 
and  the  associated  parameter  values  are  identified  simultaneously  by  minimizing  the 
fitting  residual  measured  in  the  observation  space.  The  former  requires  solving  a 
continuous  max-min  optimization  problem,  and  the  latter  requires  solving  a 
combinatorial  optimization  problem.  The  validity  and  applicability  of  the  proposed 
methodology  are  demonstrated  by  numerical  experiments.  We  have  shown  that  the 
choice  of  an  objective  function  in  model  application  impacts  the  determination  of  the 
parameter  structure  complexity. 

Additionally,  we  have  developed  a  methodology  for  optimal  observation  network 
design  for  parameter  structure  identification  in  groundwater  modeling.  The  design 
objective  is  to  find  the  minimum  cost  design  that  can  provide  sufficient  information 
for  parameter  structure  identification.  By  incorporating  the  sufficiency  requirement  as 
a  constraint  in  the  optimization  model,  the  proposed  methodology  quantitatively 
unifies  observation  network  design,  model  structure  identification  and  model 
application  reliability.  We  use  sequential  Gaussian  simulation  to  generate  a  large  set 
of  realizations  which  are  used  to  represent  the  true  parameter  field.  For  each 
realization,  we  search  for  the  minimum  cost  design  that  satisfies  the  data  sufficiency 
requirement.  After  solving  the  design  problems  for  each  of  the  realizations,  we 
analyze  the  overall  results.  We  calculate  the  average  model  application  errors  for 
different  parameter  structures,  the  probability  that  the  model  application  error  takes 
on  a  value  that  is  less  than  the  specified  accuracy  requirement,  the  minimum  number 
of  wells  required  for  the  specified  parameter  structure,  and  the  priority  of  the  well 
locations.  We  also  analyze  the  reliability  of  providing  sufficient  data  for  each  given 
cost.  The  proposed  methodology  provides  stochastic  solutions  for  the  decision¬ 
makers  as  a  reference  for  the  observation  design.  The  decision-makers  can  specify 
different  accuracy  and  reliability  requirements  and  find  the  trade-off  between  the 
requirements  and  the  design  cost.  The  proposed  methodology  can  provide  a  more 


flexible  and  robust  solution  than  the  traditional  design  methods. 
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